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Fermion Propagator in Cosmological Spaces with Constant Deceleration 
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We calculate the fermion propagator in FLRW spacetimes with constant deceleration q = e — 1, 
e = —H/H 2 for excited states. For fermions whose mass is generated by a scalar field through a 
Yukawa coupling m = gY<j>, we assume 4> oc H . We first solve for the mode functions by splitting the 
spinor into a direct product of helicity and chirality spinors. We also allow for non-vacuum states. 
We normalise the spinors using a consistent canonical quantisation and by requiring orthogonality 
of particle and anti-particle spinors. We apply our propagator to calculate the one loop effective 
action and renormalise using dimensional regularisation. Since the Hubble parameter is now treated 
dynamically, this paves the way to study the dynamical backreaction of fermions on the background 
spacetime. 
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I. INTRODUCTION 

Partly motivated by our increased capability in recent 
years to perform accurate observations on the sky, we 
have turned our attention to investigating the impact 
of quantum effects on the evolution of the universe. In 
pursuit of this goal, quantum field theory in cosmologi- 
cal spacetimes continues to be an increasingly important 
field of research. Among other things, it investigates how 
quantum fluctuations affect the background spacetime in 
perturbative quantum gravity, a process also known as 
quantum backreaction. 

A line of research deals for example with the back- 
reaction of quantum fields whose spectrum is nearly 
flat. Examples of these fields are the minimally cou- 
pled massless scalar and the graviton. Consequently, 
these fields are expected to yield a substantial backre- 
action 0, s, i, i, s, i s s a m m m m E3 ■ Ais °' 

fermions in curved spacetimes and in particular maxi- 
mally s ym metric spacetimes were investigated (see e.g.: 

El M HE 01). 

Another approach to this question is concerned with 
calculating the possible effect of the trace anomaly on 
the background spacetime. Some authors argue that this 
effect, particularly in relation with the cosmological con- 
stant problem, could be significant @, EH, S3, Hi, H, H , 
however also see [2fjj . 

An essential element to study quantum effects in 
curved spaces is the propagator. Due to its high de- 
gree of symmetry, de Sitter spacetime proves to be ideally 
suited for calculating various quantum effects. Moreover, 
results are immediately applicable to inflation, cosmolo- 
gists' favourite paradigm for a brief exponentially fast 
expansion in the early universe. Chernikov and Tagirov 
calculated the scalar propagator in de Sitter spacetime 
[27| . also see [28|, [29[. The vector propagator was con- 



structed by [3fJ, |3l| and t he g raviton propagator received 
contributions from pi. 1321 I33L l34l. l35l| . Candelas and 
Raine calculated the fermion propagator [l5l ]. 

However, de Sitter spacetime suffers from several draw- 
backs. Firstly, pure de Sitter spacetime is never realised 
in nature. Since de Sitter spacetime corresponds to a 
globally constant Hubble parameter H, this is in reality 
never attained. Secondly, de Sitter spacetime is non- 
dynamical. It is therefore inconsistent to study backre- 
action effects while at the same time assuming a con- 
stant background. We should allow the background to 
change to encompass all possible backreaction effects con- 
sistently. Therefore, it is necessary to consider propaga- 
tors in more general spacetimes. 

The scalar and graviton propa gato r in quasi de Sit- 
ter spacetime were calculated by [ll| and subsequently 
generalised to Friedmann-Lemaitre- Robertson- Walker or 
FLRW spacetimes with constant deceleration by [HI, [HI, 
[lij . In this paper we calculate the fermion propagator in 
FLRW spacetimes with constant deceleration. 

We should touch upon an important issue. In the mas- 
sive case we need two additional constraints to be satis- 
fied in order to solve for the propagator: 
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Recall that all interesting epochs in the evolution of our 
universe satisfy the first constraint. In the matter era 
we have e = 3/2, in the radiation dominated epoch we 
find e = 2. It also serves as an approximation to both 
inflation and the current dark energy dominated epoch 
when e <C 1. Note that e coincides with the slow-roll 
parameter of inflation and it is straightforwardly related 
to the somewhat more familiar deceleration parameter 
q = e — 1 . Hence (jla|) is equivalent to requiring a constant 
deceleration parameter. 

The second assumption (|lb[) is required to analytically 
derive the propagator in the massive case. In Yukawa 
theory, the mass of the fermion is generated by a scalar 
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field for which we can achieve 4> oc H in several cases to 
which we will turn shortly. 

We then calculate the one loop effective action in- 
duced by fermions using our propagator. The one loop 
backreaction arises from integrating out a free, quadratic 
fermion field and, using dimensional regularisation, this 
generates a correction to the (classical) Friedmann equa- 
tions. 

The Hubble parameter occurring in this effective ac- 
tion is now a dynamical quantity. Consequently, the 
backreaction of these fermions on the background space- 
time can thus be analysed dynamically. Also, when the 
fermions are coupled to a scalar field, this opens up the 
possibility to study the impact of the fermions on the 
evolution of scalar fields. 

This paper is outlined as follows. In section [IT] we 
review the basic theory required to study fermions in 
curved spaces and we establish our notation. Sections 
IIIII and IIVI are devoted to deriving the fermionic propa- 
gator in FLRW spacetimes in the massless and massive 
case, respectively. Finally, in section [V] we apply our 
propagator to calculate the one loop effective action. 



II. FERMIONS IN FLRW SPACETIMES 
A. The Dirac Equation 

Fermions are in general D dimensional curved space- 
times described by the action [3^| : 

s = J dW^jl [tfvrf - (v^) t'V] - ™K>} , 

(2) 

where the Dirac matrices 7^ satisfy the following anti- 
commutation relations: 

{-f,-f} = -2<r. (3) 

Variation with respect to "0 = -0^7° yields the equation 
of motion the fermion field tp satisfies: 

ry^V^VO) - mtp(x) = . (4) 

We will make use of the vierbein formalism which can be 
thought of as a transformation of the metric tensor to a 
locally flat Minkowski metric: 

9n„(x) = e%(x)el(x)r)ab , (5) 

where rj a b = diag(— 1, 1, • • • ,1) is the Minkowski met- 
ric. We use Greek letters to run over spacetime indices 
whereas we use Latin letters either for the tangent space 
indices or for spinor indices. It will be clear from the con- 
text which interpretation we are using. The Minkowski 
metric and the (flat space) Dirac matrices j a arc space- 
time independent. Spacetime and tangent space indices 
are raised or lowered by making use of the full metric or 



V(0) 




Figure 1: Qualitative sketch for the "Mexican hat" poten- 
tial fl,5\) assuming fi 2 <C \£R\- Near the minimum of the 
potential 116\) . assumption lib]) is satisfied. 

Minkowski metric, respectively. The covariant derivative 
acting on a Dirac spinor is defined as: 

V^{x) = d^(x) - T^(x) , (6) 

where the spin connection is given by: 

IV = - \e" c {d a e vd - T a ^e ad ) [ 7 C , l d ] , (7) 

such that V M 7„ = 0. Let us specialise to flat FLRW 
spacetimes in which the metric is given by g^ v — 
a 2 (jf)ri^ v . Here, a(r)) is the scale factor of the universe in 
conformal time defined by dt = a{rj)dr). In FLRW space- 
times the vierbeins are a function of conformal time only: 

e£(ri) = S^a- 1 ^). (8b) 

We thus find: 

iyV^(a) = a~ B i 1 (r))ij b d b (a 2 ^ 1 fa)V(aO) ■ (9) 

This useful identity relates covariant derivatives to par- 
tial derivatives in the tangent space. 

Since we work in the chirality representation, the flat 
space Dirac matrices are given by: 

/=(::)• m 

where I is the 2x2 identity matrix and: 

and where the o % denote the Pauli matrices. 

In Yukawa theory the mass of the fermion appearing in 
equation @ is time dependent and generated by a scalar 
field. Let us thus briefly consider the complex system 



3 



of a fermion whose mass is generated by a homogeneous 
scalar field <j>(r]), both of which arc coupled to gravity. 
The action is: 



d D x, 



1 



+ J dW^|- [i'l^^ - (V^) W] I -(12) 

Here g v the Yukawa coupling constant. By comparing 
equations (fT2"|) and @ we can thus identify the mass of 
the fermion as: 



m(r)) = g^cj>(r]) 



(13) 



Furthermore, the (classical) evolution of the homoge- 
neous background field 4>(rj) is governed by: 



□0 - d^V(cb, R) = , 
where V((f>, R) is the potential: 



V^R) 



^12 



(14) 



(15) 



Here, £ denotes some non-minimal coupling to gravity 
and R is the Ricci scalar. In FLRW spacetimes in D = 4 
we have R = 6(H+2H 2 ) = 6(2-e)H 2 , where the Hubble 
parameter is given by H = a/ a. Dots denote derivatives 
with respect to cosmic time and dashes correspond to 
conformal time derivatives. Assuming [i 2 <C £R, the so- 
lution of (fill) is: 



^/£[e(3-2e)-6£(2-e)]. (16) 



Clearly, <ft oc H and assumption (|lbj) is satisfied. More- 
over, observe we need £ < e(3 — 2e)/{6(2 — e)}. The 
potential is depicted in figure fTJ Given a potential of 
the form (|15p . we identify some interpretations leading 
to this scenario. In the early universe the Higgs field 
satisfies pi 2 <C and hence (fTB)) . A quintessence field 
has fi 2 — and (fT!))) is also satisfied in the late uni- 
verse. Finally, we can infer that a massless minimally 
coupled scalar field is also allowed in the regime where 
< e < 3/2 or e > 2. 

We work in the approximation where e = 0, which is 
not equivalent to the standard slow-roll approximation of 
inflation, where one introduces the slow-roll parameters 
e and rj. These parameters are allowed to vary slowly in 
time such that e and fj are higher order in the slow-roll 
parameters and can be neglected when one is interested 
in the leading order behaviour only. When e is allowed 
to vary in time, equation (|16|) generalises to <p — aH, 
where: 



e(3 - 2e) + 



H 



6£(2 



0(d). (17) 



The leading order slow-roll approximation to this equa- 
tion corresponds to: 



^[3e-6£(2-e)] 



(18) 



Note the equation above depends on e, whereas it does 
not depend on rj. 



B. The Propagator 

Depending on which pole prescription one uses in 
the spirit of the Schwinger-Keldysh formalism (see e.g. 
[I?], [38|, 39]), one can obtain different propagators such as 
the anti-time ordered, time ordered or Wightman prop- 
agators. We are primarily interested in the time ordered 
or Feynman fermion propagator, defined by: 

iSf(x,x) = (fi|T{^ (a;)&(x)}|n} 

= 9(r) - fj)(n\i> a (x)$ b (x)\n) (19) 

-e^- v )(n\^ b (i)Mx)\n) . 

The symbols a and b are the spinor indices and |f2) de- 
notes the state of the system. Also note the minus sign in 
front of the second ^-function. The fermion propagator 
at tree level satisfies: 

[H^l - m] iS£ b (x, x) = iS D (x - x)I ab , (20) 

and a likewise equation when x and i are interchanged. 

C. Properties of FLRW spacetimes 

We already discussed that FLRW spacetimes can be 
characterised by the scale factor a as a function of either 
cosmic or conformal time. There are however a few more 
important properties of this spacetime that are worth 
mentioning here. 

An important relation that we will use throughout this 
manuscript is: 



H( V )(1 - e) 



(21) 



This equation is equivalent to assumption pap , see [12j, 
Il3j . Let us define some relevant geometrical functions: 



y ++ (x,x) 



Ax+ + (x, x) 



'/'/ 



(22a) 



y+-(x,x) = —(- (77 — 77 + ief + \\x - x\\ 2 ) (22b) 
r\r\ V / 



y-+ 



(x,x) = -4f - ( rj - fj - ief + \\x - x\\ 2 ) (22c) 
7777 V / 

y—(x,x) = -^(-(\ V -f ] \ + is) 2 + \\x-x'\\ 2 ), (22d) 
7777 V / 
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which all vanish near the lightcone. Here, e > in (|22[) 
refers to the Feynman or time ordered pole prescription 1 . 
In de Sitter spacetime, i.e.: when e = 0, the de Sitter 
invariant function y + +(x, x) is related to the geodesic 
length l(x,x) as y ++ \ E=0 = 4sin 2 (iJ//2). 



III. FLRW FERMION PROPAGATOR: 
THE MASSLESS CASE 

Massless fermions in FLRW spacetimes are not diffi- 
cult to deal with. By making use of equation Q the 
propagator for massless (conformal) fermions should sat- 
isfy: 



Let us recall the following identity: 
1 4tt d/2 



Axl- + 2 (x,x) r(£>/2-l) 



i5 D {x - x) 



(23) 



(24) 



We can thus immediately infer the solution of equation 
(ESI): 



iS c (x, x) = (aa) 



■2ji r(D/2-l] 
4tt°/ 2 



ij b d b 



(25) 

Here, a = a(fj). This completes the calculation for the 
propagator of massless fermions in FLRW spacetimes in 
D dimensions. This propagator is valid in any FLRW 
spacetime and assumption (fla|) can be relaxed. Because 
massless fermions are conformal in any dimension, their 
propagator is much easier to calculate than the massless 
scalar propagator [ll|, [H, [H, [H| . 



IV. FLRW FERMION PROPAGATOR: 
THE MASSIVE CASE 

A massive fermion is not conformal and hence its prop- 
agator in general contains some complicated mass depen- 
dence 2 . We firstly solve for the fermionic mode functions. 
We generalise the approach outlined in [|(| El| to incor- 
porate FLRW spacetimes with constant e. Using these 
mode functions, we then return to position space to con- 
struct the Feynman propagator. 



Note that the e of the Feynman pole prescription is unrelated to 
the "slow roll" parameter e in FLRW spacetimes. 
Note in [3r3 | there is an erroneous statement regarding massive 
fermion propagators in curved spacetimes. It is argued on page 
87 that a fermionic propagator in any spacetime can be related to 
the scalar field propagator in that spacetime. This is not correct. 
The (spinorial) structure that arises when a covariant derivative 
acts on a spinor (01 is much more complicated than for scalar 
fields. 



A. Fermionic Mode Functions 

Let us firstly define the rescaled fermionic spinor: 

x (x) = {ri)i>{x) . (26) 

Keeping an eye on equation ^ one can easily check that 
this factor is chosen such to conveniently transform the 
covariant derivative into a partial derivative. Equation 
of motion (j4j thus reads: 



ij b dbx( x ) ~ am x( x ) = • 



(27) 



The canonical momentum associated with ip{ x ) follows 
from promoting the variation of ^ with respect to i/j{x) 
to an operator. According to the usual moves in quan- 
tum field theory, we impose anti-commutation relations 
between ip(x) and its associated canonical momentum: 

{VUx, t),a D -\t)r b {y, t)} = 6 D -\x - y)6 ab , (28) 

with the other anti-commutators vanishing. By making 
use of the rescaling in equation |26|) , appreciate that the 
anti-commutation relations above simplify: 



{xa(x,t),r b (y,t)} 

{Xa(x,t),Xb(y,f)} 

{r a (x,t),xt(y,t)} 



S D -\K-y)S ab (29a) 
(29b) 
. (29c) 



1. Chirality and Helicity Decomposition 

Therefore we expand the rescaled spinors x( x ) an< i 
x(x) in creation and annihilation operators as follows: 



x( x ) 



x( x ) 



PE«"^(M)e ta (30a) 



(27T) 



(30b) 



&l h x {h) (Kv)c 



-ik-x 



The expansion above merits a few remarks. Firstly, a^h 
and fek./i are the fermion and anti-fermion annihilation 
operators of helicity h, respectively, in the usual sense: 
ak,h|^) = = 6k, h\ty- The helicity h, i.e.: the spin in 
the direction of motion, can be either +1 or —1 in units 
of h. 

Let us for simplicity return to the D = 4 setting we 
are familiar with. We will generalise the following con- 
siderations shortly to arbitrary dimensions. 

In the equation above x^^C 4 ; 7 ?) i s a 4-spinor of helic- 
ity h. We decompose the 4-spinor x(k, rf) into a direct 
product of chirality and helicity 2-spinors: 



x(M) = £ x (ft) (M) = E 



x^hiKv) 
XrAKv) 



(31) 
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Here, XL,/i(k, rj) and XR,fc(k, 77) are left- and right-handed 
1-spinors of helicity h, respectively. Furthermore, ^ is 
the helicity 2-cigenspinor: 



(32) 



where a = (a 1 ,<j 2 ,cr 3 ) is a shorthand notation for the 
Pauli matrices. Furthermore note that h 2 = h^h = 1. 

For future convenience, let us explicitly derive the two 
vectors We write: 



k — (kxi ky,k z ) . 
and keep in mind that: 



I ^11 2 — k 2 , + k 2 + k 2 — 1 



(33) 



(34) 



It is a trivial exercise to solve for the two eigenvectors of 
h, which correctly normalised to unity read: 



1 



2(1 -k z ) 
1 

2(1 + fc,) 



h x thy 

i-L 



iky k a 

i + L 



(35a) 
(35b) 



Moreover, the helicity eigenstates are mutually orthogo- 
nal: 



0. 



(36) 



Finally, we have tacitly ignored the expansion of the 
anti-particle contribution f(k, rj) in chirality and helicity 
spinors in equation (|30p . We will return to this subtlety 
shortly. 



2. Generalisation to Higher Dimensions 

Our chirality and helicity decomposition of the spinor 
degrees of freedom in equation (|3"Tj) also works in D space- 
time dimensions. Once a particular representation of the 
gamma matrices has been found, the projection operator: 



P+ = 



1±7 Z 



(37) 



splits the 2 D/2 -spinor degrees of freedom in two equal 
contributions of definite chiralit y, i .e.: a left- and right- 
handed 2( D ~ 2)/2 -spinor (see e.g. 021 )• Here: 



7° +1 = a D"f° ■ ■ ■ 7° 1 



(38) 



where a D is fixed by requiring that P± is a proper pro- 
jection operator: 



P±=P± 



This yields: 



P+P- = 



(7 D+1 ) 2 = 1 



OLt 



exp 



i-{D-l){D + 2) 



(39) 
(40) 



In D = 4, we recognise the familiar 7 s = i7°7 1 7 2 7 3 . 
An helicity operator that commutes with the projection 
operator can be defined in Fourier space as: 



h = — 



. (P-2)(P-1) 



T(D-l) 



Alt— It 



2 ka h ---r 



(41) 



where e" 1 " ' 1d ~ 2 is the Levi-Civita tensor in D — 1 di- 
mensions. The exponential phase ensures that h is her- 
mitian. Finally, the Gamma function has been inserted 
to account for the number of different permutations of 
gamma matrices. 

Hence, X L ,fe(k,?7) and XR,/i( k , rj) are 2 (D ~ 4)/2 -spinors. 
Note the definition of the helicity operator is consistent 
with (|3"2")) . Most importantly, the chirality and helic- 
ity decomposition as proposed in equation (|31[) carries 
through for any dimension. 

When an explicit form of £/j is required, we return to 
D = 4 and use equation (|3"S"j) . We subsequently generalise 
the result of such calculations by analytically extending 
it to arbitrary D. 



3. Spinorial Normalisation Conditions 

In textbooks on quantum field theory, such as [43|, [44| , 
one expands the fermion in momentum space in spin 
eigenstates, whereas we expand in helicity eigenstates. 
This results in a different normalisation requirement on 
the spinors in momentum space and consequently the 
standard textbook results cannot be straightforwardly 
copied. The reason for expanding in helicity eigenstates 
is that in curved spacetimcs helicity is more convenient 
to work with rather than spin. 

Spin is a conserved quantity, i.e.: an appropriate quan- 
tum number, in static backgrounds, whereas helicity is 
conserved in time dependent but spatially homogeneous 
backgrounds (4f|. One can easily show that the helicity 
operator constructed above in (|4ip commutes with the 
kinetic operator in Wigner space for spatially homoge- 
neous correlators. 

We normalise the spinors using two conditions: a con- 
sistent canonical quantisation and orthogonality of par- 
ticle and anti-particle spinors. 

We thus impose the usual anti-commutation relations 
between creation and annihilation operators: 

{a k ,fc,4',fc'} = (27r) D - 1 ( 5 D - 1 (k-k')<5 /l ,^ (42a) 

ik,h,bl,h'} = (27r) D - 1 5 D - 1 (k-k')<S, l ,v, (42b) 

and all other anti-commutators vanish identically. Now 
(|42]) is consistent with ((29)) if and only if: 



E ( k ' ^l (h) ( k ' v) + "i h) (-k, vH {h) (-k, n) = s ab . 

h 

(43) 

Note that this is a condition on 2 { - D / 2 ^~ 1 {2 D / 2 + l) matrix 
elements in spinor space, i.e.: for each a and b the above 
equality has to be satisfied. 
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Secondly, we require that the particle and anti-particle 
spinors are mutually orthogonal: 

X W (k,rj)^ h '\k,r]) = = D {h ' ) (k,r 1 )x {h) {k,?]). (44) 

This restricts 4 matrix elements in helicity space, i.e.: for 
each h,h' = ± the above equality needs to be obeyed. 

Pcskin and Schroeder [43l | impose normalisation con- 
ditions on both particle and anti-particle spinors in spin 
space. The corresponding conditions in helicity space are 
only partially satisfied. Both particle and anti-particle 
spinors of different helicity are trivially orthogonal by 
construction ()36|) . However, we do not relate the am- 
plitude of equal helicity states for particles and anti- 
particles to the mass of the spinor because helicity eigen- 
states for a particle at rest cannot be constructed. 



4- Particle Mode Functions 

We will now solve the Dirac equation. Let us insert 
equation ([50)1 into the Dirac equation (|27[) . Making use 
of our chirality and helicity decomposition (|31 [) trans- 
forms the Dirac equation to a first order coupled system 
of differential equations: 

+ hkx^,h - amxaM = (45a) 
*Xr,^ - hk Xn,h - amxL.h = 0. (45b) 

The arguments of the functions in the equation above 
have been omitted for notational convenience. The left- 
arid, right-handed spinors however are now a function of 
the magnitude of the Fourier mode k only, because the 
differential equation is invariant under k — > —k. 

Standard model fcrmions are chiral, in the sense that 
m R =/= m L . For simplicity, we take one mass m = m R = 
m L only. The following derivation can easily be gener- 
alised when standard model fcrmions are considered by 
making use of the P± projectors in equation (1371) . 

We now return to the expansion in chirality and helic- 
ity spinors of the anti-particle contribution ^(k, rf), which 
we did not need so far. Since the Dirac equation is a lin- 
ear differential equation, the solution for i/(k, rf) cannot 
contain any new degrees of freedom once we completely 
solved for x(k, rf). In other words: the form of v(k, rf) is 
dictated once we have completely specified x(k, rf). We 
expand i/(k, 77) slightly differently: 



KM) = $>W(M) = £ 



v^ h {k,rf) 
fL,fe(k, rf) 



(46) 



Note we flipped the position of the left- and right-handed 
spinors compared to (|3Tj) . Consequently, the resulting 
equations of motion for both ^ R ;,.(k, rf) and zy L! /j(k, 77) are 
identical to ([25]). 

Alternatively, we could have expanded differently and 
flipped helicity, i.e.: send h — ► —h in equation (|3T|) . Be- 
cause the mass enforces mixing between different chirality 



states, labelling these states differently has no physical 
relevance. 

We return to equation of motion (j45|) and define: 

u±h(k,v) = ■ (47) 

This transforms equation ([45]) to: 

iu' +h + hkii-h — amu+h = 
iu'-h + hku+h + amu-h = . 



(48a) 
(48b) 



We exploit assumption (fla[) , which is equivalent to ([21 
Inserting this expression into the equation above yields 



hku-y 



u+h 



hk 



u +h 



H(l-e) V 
m U—h 



= 
= 0. 



(49a) 
(49b) 



Recall that the mass of a fcrmion in Yukawa theory is 
time dependent in FLRW spacctimcs. There does not 
exist a rescaling of the functions u±h that removes the 
mass dependence in the equation of motion above. Even 
when one assumes that the mass m is time independent, 
one cannot find such a rescaling. It is particularly simple 
to choose m/H constant as following from equation {T( 
and for convenience we define: 



H(l-e) 



(50) 



We can proceed and decouple the two linear differential 
equations to find: 



L ±h 



+ U 2 + 



c 2 ±< 
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u± h = . 



The solutions are given by: 



(51) 



(52) 



where a± k are two normalisation constants that still need 

to be determined and Hv± is the Hankel function of the 
first kind of order: 



This implies: 



v+ + V- = 1 

xA = v- 



(53) 



(54a) 
(54b) 



The Hankel function of the second kind also solves equa- 
tion ([51]) . In the infinite asymptotic past, equivalent to 
the deep UV, the fcrmions become effectively masslcss for 
e < 1 and the distinction between u + h and U-h vanishes 
as a consequence. Therefore we require that the vacuum 
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mode functions in this regime equal the standard confor- 
mal vacuum solutions: 



lim u+hi-kij) = ^7=e lkv 



lim U-h(—hr)) 



// — ' — DC 



-h 



,-ikri 



(55a) 
(55b) 



This excludes the Hankel function of the second kind for 
the moment from contributing. Away from the vacuum, 
we should of course also incorporate the second solution 
to (|5ip to allow for mode mixing. Note that the UV 
limits above yield the familiar UV behaviour for the left- 
and right-handed spinors (e.g.: [H[): 



lim XL.feCM) = e ikv 

rj—^—oo Z 



lim Xr,/»(*j^) 

r}— * — oo 



l + h 



— ikrj 



(56a) 
(56b) 



When solutions l|52p are substituted in (|49ap one finds: 
a_ k = iha +k e . (57) 



We are free to fix the remaining coefficient arj_ k by the 
following condition: 



^2 \ u ±h\ 2 = \ u +h\ 2 + \u-hl 2 = 1 



(58) 



We postpone the motivation for this requirement to sec- 
tion IIV A 61 For the moment, let us simply normalise 
accordingly for convenience. We thus find: 
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(59) 



where we have chosen the unobservable phase factor in 
accordance with the UV limit (|55|) . We use some well- 
known identities involving Hankel functions which we list 
for convenience in appendix lAl 

Concluding, we find that the vacuum solutions to equa- 
tion (|49p neglecting derivatives of e and treating the quo- 
tient m/H as time independent are given by: 



u +h (-k V ) = e if(" + +V2) (60a) 
(-k V ) = -he^- +1 'V ^-^p (-lev) -(60b) 



U-h 



As we send e — > and approach de Sitter space, our solu- 
tion is in accordance with 41| , as can be seen from (|53p . 
In this solution we have tacitly assumed that e < 1 such 
that T] < 0, equivalent to an accelerated expansion. If 
however e > 1 or 77 > 0, which corresponds to a decelerat- 
ing universe, one should simply use relations (|Ald|[ATcl) 
in appendix The case e = 1 requires special attention, 
see e.g. [13] for the analogous scalar propagator case. 

Let us now turn our attention to finding the solution 
for the second mode function corresponding to the ex- 
cited state. The procedure is completely analogous: we 



again let 77 — > — 00 and match to the appropriate UV 
expansion involving the second solution Hi,±{—kri). We 
normalise the solution according to ([55]) and choose the 
phase factor in accordance with the corresponding UV 
limit. We can thus construct the total solution by linear 
superposition: 



u + h{-krj) 
u-h{-krj) 



cH( V ,C)+dH( V ,0 
-hcH*(r),0 + hdH*(r),<;). 



(61a) 
(61b) 



where: 



K.++1/2) J w( i 



H( V ,0 = e -^+^^^Hll\-k V ).(62b) 

Note that we have omitted all k dependence in ([6"Tj) for 
notational convenience. Note finally that the two funda- 
mental solutions (|62p in u+h(~ kif) are not complex con- 
jugates of each other, unlike for example for the scalar 
field case. However, the second solution U-h(—krf) can 
be obtained from complex conjugation of the fundamen- 
tal solutions comprising u+h(—krj), as can clearly be seen 
in equation (|6ip above. 

For the moment, we leave the normalisation constants 
c and d undetermined. Appreciate we are nevertheless 
free to normalise the fundamental solutions to unity. 
Of course we will determine the normalisation constants 
shortly by applying our analysis performed in subsection 
IIV A 31 But first, we need to solve for the anti-particle 
mode functions. 



5. Anti-particle Mode Functions 

Let us stress again that the other spinor u^ h > (k, 77) can- 
not contain any new degrees of freedom because Dirac's 
equation is linear. We only have to determine the con- 
nection with x (k, 77). When we expand according to 
(|46p the coupled linear differential equations ^ Ri /j(k, rj) 
and i/ L] /i(k, ??) obey are identical to (|4"5)) . Therefore, we 
can analogously define: 



v±h{k,ri) 



V2 



(63) 



We are rewarded for our extensive discussion of chirality 
and hclicity. The solutions of v±h{—krf) are identical to 
([6"Tj) . where we should only denote the coefficients in this 
solution differently: 

v +h {-krj) = fH(r),()+gH(ri,0 (64a) 

u-h(-kv) = -hfH*(r ] ,0 + hgH*(r ] ,0- (64b) 

The functions appearing in this equation are given by 
([62]) . We need to determine how the coefficients c and d 
of u±f l (—krf) relate to / and g of v±i l (—kr)). Moreover, 
we need to find the correct normalisation condition for c 
and d. 
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6. Normalisation of the Mode Functions 



Inspired by our discussion in IIV A 3[ we impose two 
normalisation conditions on the spinors: a consistent 
canonical quantisation (|43[) and orthogonality between 
particle and anti-particle spinors (|44p . 

After quite some work the, naively, 10 constraint equa- 
tions in D = 4 in (|43]) yield the following 3 conditions on 
our coefficients: 



|c| 2 
|c| 2 



Ml 2 
Ml 2 



l/l 2 
l/l 2 

cd* 



■\9\' 

far 



(65a) 
(65b) 
(65c) 



The orthogonality condition (|44| yields one more condi- 
tion: 



The solution is: 



l/l 
\fj\ 



cf* + dg* 



\d\ 

ki 



o. 



±7T. 



(65d) 



(66a) 
(66b) 
(66c) 



The last line relates the phases of our normalisation con- 
stants, e.g.: (j>f denotes the phase of /. Moreover, from 
(l65l) we can derive: 



l/P 



\9? 



(66d) 
(66e) 



where the second line is a consequence of the first. Note 
that this condition also follows from charge conservation 
as argued in 



4ll | . Changing the normalisation in l(5"8 



results in a change of the right-hand side of equations 
(|65a| . (|66djl and (|66c[) . Normalising the fundamental so- 
lutions to unity as in (|58[) is particularly convenient as it 
allows for a particle interpretation of \d\ 2 and |/| 2 . 

The normalisation constants provide us with three 
physical degrees of freedom. Of course \c\ is physical 
which determines |d| through (|66d[) . Hence, 4>d is phys- 
ical just as in the scalar field case. Since <f> c is a phase 
that can be removed without physical consequences, the 
phase relation (|66cp determines 4>f in terms of the third 
physical phase 4>g ■ Note finally that one of the two phase 
relations (|65c|) and (|65d|) is redundant in this part of the 
analysis. 

A final remark is in order. To derive (|65p we needed the 
explicit form of ^ in equation (|35p . which is only valid 
in D = 4. However, we analytically continue this result 
to arbitrary dimension and one naturally finds (|65p . 



7. Summary 

Using the solutions for u±i l (k,rj) and i J ±h(k,rj) given 
in (fBTj) and respectively, we recall relations (|4"7| and 



(|63j) to find: 

Xn,h(k,v) = 
and moreover: 
^,h{k,ri) = 
v R ,h{k,ri) = 



1 

71 
i 

71 



i 

71 
i 

71 



c{H-hH*} + d{H + hH*} (67a) 
c{H + hH*} + d{H - hH*} , (67b) 

J{H-hH*}+g{H + hH* }] (67c) 
f{H + hH*} + g{H - hH*}] , (67d) 



where the fundamental solution are given in equation 
(|62|) . The normalisation conditions (|66|) determine / and 
g in terms of c and <i Moreover, we have |c| 2 + |d| 2 = 1. 



B. The Feynman Propagator 

Having discussed the fcrmionic mode functions in some 
detail, we can turn our attention to solving the Feynman 
propagator for fermions. We insert the rescaled spinors, 
expanded in terms of creation and annihilation opera- 
tors, into the formal definition of the Feynman propaga- 
tor (fT9]l . which yields: 



(68) 



f D-l D— 1 

9(7) — fj)a 2 2 



(2tt 



£5-1 D-l 



-^(7/ — i])a 2 a~ 



ik- (x— x) 



where a = 0(77) and a = a(fj). The reader can easily see 
that we have to make use of the solutions (f6"7| above to 
calculate the 16 matrix elements of the propagator. 



1. Outline of the Calculation 

Let us outline this rather cumbersome calculation in 
seven steps. 

1. We recall the form of e.g. the massless fcrmionic 
propagator (|25p and realise we have to extract an opera- 
tor acting on x, rather than on x. 

The reader can easily see that the 16 matrix elements 
of the propagator contain products of two Hankel func- 
tions. Some of these products involve Hankel functions at 
the same order, but others encompass Hankel functions 
where the order differs. Since we can only easily integrate 
products of Hankel functions of equal order we need to 
transform precisely the latter products by making use of 
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identities (|A1[) . e.g. 



where: 



H^(-k V )Ht\-kfj) 



e ±i7r(i/_ — 1) 



d 
(1// 



V 



(69) 



where + or — applies when m equals 1 or 2, respectively. 
Clearly, the "extracted operator" is a function of 77. 

2. We rewrite the Hankel functions in terms of Mac- 
Donald functions. By making use of 0, [4?J we find: 



2i 



2i 



Hf>(z) = -e^ v K v (iz), 
ir 



(70a) 
(70b) 



where K v (z) is a MacDonald function. 

3. The emerging matrix structure of equation (|68[) 
is most easily uncovered by evaluating its 2 x 2-matrix 
constituents successively and writing the result in terms 
of Pauli matrices. Recalling equations (|3"3"j) and (f3~4")l we 
derive: 



kyOy 



k z cr. 



kj(J 3 



±i<r 3 dj 



(71) 



where the — or + sign applies depending on whether one 
deals with the 9(rj — fj) or 9(fj — 77) contribution in (f6"5| . 
respectively. The last step is possible because only the 
exponent depends on x. Note finally that the index j in 
equation (|7ip is supposed to run over (flat) spatial indices 
exclusively. 

4. The k independent contribution to the propagator 
can thus be pulled out of the Fourier integral. These 
matrices simplify considerably when we consider the fol- 
lowing identities: 



(«7 dt + am 
1 

~~ 2 



1± 7 ° 



±id n ± ia J dj 



(72) 



id v — ia 3 dj ± am ±i9 l; =F ia 3 dj 



ia 3 dj ± am 



am 



Clearly, we recognise the complex conjugate of the Dirac 
operator and the (l±7°)/2 structure we are familiar with 
from the de Sitter propagator. 

5. Let us present an intermediate result that already 
captures much of the final structure of our propagator. 
For brevity, let us only consider the contribution to the 
propagator arising from the first ^-function in the prop- 
agator, i.e.: when 77 > fj. This part of the Feynman 
propagator reads: 



D-l D-l 

a 2 a 2 



(i"/ h db + am) 



(73) 



+ 



1-7 r d^k 



(2 



—-K + (k,r),f))j\, 



K-(k,ri,fj) = \c\ 2 K v _(ikri)K v _{-ikfj) (74a) 
—icd* K v (ikrj)K v _ (ikfj) 
—ic*dK l ,_ (—ikri)K u _ (—ikfj) 
-\d\ 2 K v _(-ikrj)K v _ (ikfj), 

and where: 

K+(k,r],fj) = |c| 2 K u+ {ikr])K v+ (-ikfj) (74b) 
+icd* K v+ (ikrf)K v+ (ikfj) 
+ic*dK v+ (—ikri)K„ + (—ikfj) 
-\d\ 2 K v+ (-ikr))K u+ (ikfj) . 

The contribution to the propagator when 77 < fj is given 
by the same expression, where we only need to replace the 
normalisation constants c and d by / and g, respectively. 

The complicated structure of the propagator in spinor 
space has dramatically simplified by extracting the ap- 
propriate operators in position space, which could al- 
ready be expected from the propagator in de Sitter space- 
time. 

6. We need to perform the Fourier integrals over the 
functions K_(k,mfj) and K + (k,mfj) next. We use: 



d^k 



(27T) 



-e^f(k) 



(75) 



(4ir) 2 Jo 



' J d-3 (kx) 
dkk D - 2 2 D 3 f(k) . 



which is valid for any function f(k) that depends solely 
on k = ||k||. Here, J^(kx) is a Bessel function of the first 
kind. The reader can easily verify that all the Fourier 
integrals we need to perform are of the following form: 



dkk^J^c^K^a^K^lBk), (76) 



where the arguments of the MacDonald functions are 
purely imaginary: 



a 



±irj 
±ifj . 



(77a) 
(77b) 



We perform these integrals in appendix [B] by making use 
of [4g|. The result of this integral can be expressed in 
terms of the Gauss' hypcrgcomctric function 2-Fi- The 
^-functions in the propagator in combination with the 
arguments of the hypergeometric functions allow us to 
rewrite these arguments in terms of the y-functions (|22j) 
we introduced in section HTCl 

A well-known fact in Minkowski or de Sitter spacetime 
extends to our analysis in FLRW spacetimes: the small 
real contribution needed to make the Fourier integrals 
convergent determines the e pole prescription of the var- 
ious contributions to the propagator. 
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7. The final result for the constant e and constant £ FLRW fermion propagator reads: 

(ar] afj) 



where: 



iS±(x, x) 



iSf{x,x) =a(i7% + m) 

r(f ±z C )r(^ T <) 
(^) fl/ T(f) 



'ao 



1+7° 1-7° 
iS+(a;,5) — \-iS-{x,x) — - — 



,'D D-2 D y ++ (x;x) 

\c\ 2 Fi ( — ±i(, — — t<; y;! ^ 



D-2 



Ml 2 2 F 1 (-± l c 1 ^—TiC;^ 1 - 



y — (x;x) 



(78) 



(79) 



where the y- functions are given in equation Q22|) and £ in (|50p . Moreover, to write the propagator in terms of the 
y-functions we introduced the following notation: 



x = (fj, x) = (-??,x) 
I 



(80) 



One can think of x as an antipodal coordinate (see e.g.: 
[IS IIS B El, [13 )■ Note that we extract an operator in 
([75)) that does not depend on x. However, there is noth- 
ing special about the x leg of the propagator. Keeping 
remark 1 above in mind, we know that we could have 
equally well written the propagator (|78[) in terms of an 
operator acting on x by replacing: 

a (H M V M + m) -> a (i7 M V M + rrij . (81) 

This propagator would then have satisfied the Dirac 
equation (|20p in terms of x. 

As a first check of this result, note that if we send 
to — > and d — ► 0, the massive fermionic FLRW propa- 
gator ([75]) correctly yields the massless fermionic FLRW 
propagator earlier derived in equation (|25p . Again, in or- 



der to derive this result we needed the explicit form of £h 
given by (|35|) . The operator in (|78|) has thus been con- 
structed in D = 4. However, we can analytically continue 
again to arbitrary dimensions. 

Finally, let us study the IR behaviour of our propa- 
gator. Note that the coincidence propagator (in Fourier 
space) is IR finite due to the Pauli exclusion principle. 
In position space, note that the IR limit of all the y- 
functions in equation (|22[) is: 

y(x;x)^±cx>, (82) 

where the plus or minus sign applies for infinite space- 
like or timclikc separation, respectively. If we recall the 
following well-known relation: 



T(c)T(b-a) f , t ._ , 1 A , T(c)T(a-b) M „ 1 



2fl(a ' iiC;z) = #4 (1 " Zr ° rt^c-bia-b + l;— J+ r ; a y (c _ 6 S (l-*)- 6 ^c-aib-a+V,.^, 

(83) 

we can easily extract the leading order IR behaviour of the propagator in equation (|T9"|) : 

iS±(Iii) . r (w )r(1±2iC) (i<f rstttesy'^ T «. e<-^ ffcM) v ( ^ ±i£) (84) 



(47r) D/2 r(l±<) [ V 4 / V 4 



Note that for the spacclikc IR, all y-functions behave identically and the result above can be simplified further. 
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For the timclikc IR, the e-prcscription in equation (|22[) 
dictates how to take the powers of the y-functions. Of 
course, one has to insert equation (|84|) above into (|78|) to 
obtain the full propagator in the IR limit. 

The fermionic propagator (f78|) is analogous to the 
scalar and graviton propagator in universes with constant 
deceleration parameter q = e — 1 derived in [12, EH, EI| • 
Unlike in the case of scalars and gravitons, where quite 
generically the Bunch-Davies vacuum is IR divergent, the 
fermionic propagator is IR finite due to the Pauli exclu- 
sion principle which forbids an accumulation of fermions 
in the deep IR. The IR behaviour for the scalar and gravi- 
ton are important, as it can comprise secular effects (that 
gro w as a power in time) in for example T4ff($) (see: 
[5l|) and (see: [HI, [EH). Mathematically speaking, 

Pauli blocking is enforced by the imaginary i in equation 
(|5ip . which ultimately carries through in the complex in- 
dices of the hypergeometric functions in the propagator 

GS). 



C. Fermionic de Sitter Propagator 

Having derived the fermionic FLRW propagator, let 
us verify whether it correctly reduces to the fermionic 
de Sitter propagator, which is a well-known result in the 
literature. The fermion propagator in de Sitter spacetime 
reads [lEGl]: 



iSp(x, x)=a (i7 M V M + m) 



(85) 



H L 



■ qa i ~ ^+r , . oA/ ^ 1 - r 

ib + (x, x) — 1- ib_(x, x) — - — 



where the superscript A denotes dc Sitter spacetime and 
where: 



iS±(x, x) 



i r(f -i Ti 2)r(f ±zf) 



' D i 



(4tt)" /2 r 

„ I D m D m D v 

X2^i it i—, — ± i— ; — ; l - - 

V 2 T if' 2 H' 2 ' 4 



(86) 



Here, y is shorthand for y ++ {x,x). If we send e — > 
and set c = 1 and d = 0, equation ([78]) above correctly 
reproduces the de Sitter result. Moreover, using: 



D 



D D y 
1 — • — • 1 - - 

' 2 ' 2 ' 4 



4/ 



(87) 



note that the masslcss limit of (j85|) indeed corresponds to 
(|25|) in dc Sitter space. We conclude that the fermionic 
propagator in FLRW spacetimes correctly reduces to the 
known cases in existing literature. 



to the effective action. This is important to study the 
impact of fermions on the evolution of the background 
spacetime and of scalar fields coupled to these fermions. 
We renormalise using the minimal subtraction dimen- 
sional renormalisation technique. 



A. Evaluating the One Loop Backreaction 

The one loop effective action formally reads: 
ri = -iTrlog[V=5(i7"V M -m)] 

/m 
dmTr[ y /=giS#'(x,x)] . 



Here, the trace is both over spatial and spinorial indices. 
We have been able to evaluate the logarithm at the ex- 
pense of losing all mass independent contributions to the 
effective action. 

The effective action for a massless fermion is com- 
pletely fixed by the trace anomaly which this field is 
known to exhibit (see for example [25|, [2(| [36| ) . The trace 
anomaly can be generated from a finite non-local effec- 
tive action (see e.g.: [25[). It can alternatively begener- 
ated from an infinite but local effective action [26|. We 
will henceforth neglect all contributions from the trace 
anomaly to the effective action in this paper because it 
has already been extensively discussed in the literature. 

We treat the covariant derivative according to © and 
we can easily expand the hypergeometric functions in 
equation (|79j) at coincidence. In dimensional regularisa- 
tion, the D dependent powers in this expansion do not 

contribute. These arc generated for the y++ and y 

contributions to ([7^]) by the first term in: 

X2F1 (c — a, c — 6; c — a — b + 1; 1 — z) 
T(c)T(c- a-b) 



T(c-a)T{c-b) 



1F1 (a, b; a + b — c+ 1; 1 — z) . 



Moreover, when we trace over spacetime indices note 
there is no contribution at all from the derivative term 
hitting one of the hypergeometric functions: 



lim > 7 b <9t y T t + (x, x) = — lim > 717' 

x — *x — * x — >x — ' 



l y++{x 1 x) 



n=0 



n 



0. 



(90) 



Here, n is the integer valued coefficient of the Taylor ex- 
pansion of a hypergeometric function. A similar identity 
holds for y Likewise, for y_| or y |_, we have e.g.: 



ONE LOOP EFFECTIVE ACTION 



lim "/ b d b y+-(x,x) = 0. 



(91) 



As a simple application for the propagator derived 
in section IIVB1 we calculate the one loop contribution 



Note that both and y h depend on the antipodal 

coordinate x in equation (|78[) . At coincidence, we also 
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use: 



The result is: 



(aa) 2 (r/fj) 2 am =F 



D-2i 
2 77 



(92) 



±i(H\l-e\y- 2 H(l-e){^ T 



dW ^^P jr (1 - 1 ) (| C | 2 - |<f ) 



- e) fc 



r( f + iC r f-ic 



(1 - ic) r (1 + »c) 



+ die" 



dm ■ 



iC r 



r(f) 



r 



(93) 



Let us compare this result with the known calculations 
in the literature in de Sitter spacetime (e — ► 0) and in 
the vacuum (\d\ = 0). Note the two omissions in the 
Candelas and Raine effective action [l5[- When tracing 
over spinor indices, we have: 



Tr 



1± 7 ° 



(94) 



and this trace does not, as is apparent from [15j, equal 
D/2 because a spinor in D spacetime dimensions has 
2"2 degrees of freedom (see e.g.: HH). Indeed, this has 
already been noted by [53|, [5J]. Moreover, the effective 
action in [l5[ misses a factor of 2 stemming from the 
separate particle and anti-particle contributions, which 
has been corrected for in e.g. [l9l ]. 

The reader can easily verify that (|93l) , taking these er- 
rors into account agrees with [HI, Il9j in de Sitter space- 
time and in the vacuum. 



B. Dimensional Regularisation 



we we neglected 0{(D — 4) 2 } contributions. Finally, we 
recall the familiar result: 



r 1 



D 



D-4 



7e - 1 + 0(D - 4) , (97) 



where 7 E is the Euler-Mascheroni constant. We thus 
arrive at: 



r, = 



d ^[(N 2 - 

\2 v ' 4 



2\ ,,D-4 



D-4 



■ 7a - 1 + In 



H 2 (l-e) 2 



(98) 



+ / dm(m 3 +mH 2 (l-e) 2 )\?P(2-iC)+ip(2+i() 

r m 

+2H 3 (1 - e) 3 Im (cd*) / dm T(2 + <)r(2 - iQ 



In the spirit of dimensional regularisation, we now ex- 
pand around D = 4. We can use the familiar result: 



T(x + ye) = T(x) [1 + ye x/j(x)} + 0(e 2 ) . 



(95) 



where e <C 1 and where ip(x) is the digamma function 
defined by ip{x) = d\ogT(x)/dx. Anticipating the form 
of the counterterms we will add shortly, we introduce a 
scale \i by: 



(H\l-e\r = (H(l-e))Y 



H\l-e\ 



(96) 



(ff(l- £ ))V 



D-4, fH 2 (l-e) 2 
1 H In ' 



Here, Im(cd*) = |c||d| sin(t/) c — <fid) denotes the imaginary 
part of cd* . 

The one loop backreaction contains a divergence when 
D = 4. However, appreciate that the contribution to the 
effective action multiplying the mixed coefficients is finite 
in D = 4. 



C. Tree level Friedmann Equations 

Anticipating the renormalisation procedure in the next 
section, we need the tree level, i.e.: classical, Friedmann 
equations of motion responsible for driving the expansion 
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of the universe. The system of interest is thus given by: promote this constant to a dynamical quantity: 



S 



- I 

K 



d°Xy 

d D x, 



-g{R-(D-2)A} 
1 



-fJ 



(99) 

gVBrtpdvtp-Vfa)} 
g^d^d^ - 1/(0) - g^Wt 
[^V^ - (V^) W] 



Here, k = 167tGn represents the rescaled Newton con- 
stant and A denotes the cosmological constant. The 
reader can easily recognise the usual Einstein-Hilbert ac- 
tion on the first line. We introduce a new scalar field 
ip, responsible for the dynamics of the universe. By an 
appropriate choice of the potential V(<p), a scalar field 
can mimic any mixture of fluids relevant for the evolu- 
tion of our universe (see e.g. |55|). The last two lines in 
equation (f9"9")l above are identical to (fT!?|) and contain the 
fermion field and the scalar field <j> generating the mass of 
the fermion through a Yukawa coupling, as in (fT3j) . We 
write: 



ip(x) = <p (t) + 5<p(x) 
</>(x) = (j> Q {t) + 6(/)(x) , 



(100a) 
(100b) 



where we assume the background fields (po(t) and 4>o(t) 
to be homogeneous. Moreover, we assume that p v ^> p^, 
such that <p(x) drives the dynamics of the universe as 
stated above. The classical Fricdmann equations of mo- 
tion are: 

... 1 1 , 



H 2 - 



1 



D - l 



H 



D - 1 



-H 2 



-A- 
A 

2 + 



V(<po)) 



(U-l)(D-2) 
«(i0g-V^o)) 



2(D-2) 



= (101a) 
= 0,(101b) 



and moreover, the classical scalar field equation of motion 
reads: 



dV 

<p + (D-l)H<p + —(tp ) 



= 0. 



We can thus derive the following identities: 

r dv ( \ 

>K — (<A)) 



y/2(D - 2)eH 



dip 



= -y/2(D-2)e(D-l- 



'Ktpo 



a 2 v 

dip 2 



-^j2{D-2)eeH 2 
2(D-1- e)eH 2 . 



(101c) 

(102a) 
s)H\l02b) 
(102c) 
(102d) 



Here, wc have used assumption (ITal 



D. Renormalisation 



In order to derive equation (|98|) . we have assumed a 
constant deceleration e, as in equation (fTa|) . We now 



5(t). 



(103) 



We will motivate this step shortly. The effective action 
(f9"5| now contains divergences in both H (t) and e(t) that 
we ought to cancel by an appropriate counterterm action. 
It is not possible to identify local covariant counterterms 
of curvature invariants only that remove the singularities 
in both of these quantities at the level of the effective 
action. We can use the tree level equations of motion de- 
rived above in (|102[) to rcnormalise the effective action, 
which is first order in h. By making use of an infinitesimal 
field redefinition, one can show the following: substitut- 
ing the zeroth order equation of motion into the first or- 
der contribution to the action yields, up to one loop, the 
same equation of motion as directly varying the action. 
We thus add the following counterterm action: 



T ct = J d D x^/~ r g{ci(j) 2 R + c 2 (, 



,d 2 V 

d<p 2 



(¥>o) (104) 



Kffg^dftipdvip + c 4 0' 



C3 



The form of the counterterm action is unique, which can 
be shown by making use of dimensional analysis and by 
requiring that only two time derivatives can act on tpo (t) 
or a(t) for stability. These conditions limit us to the 
counterterms above and possibly <p 2 ipnip. However, one 
can easily verify that the latter term does not have the 
correct form required to cancel the divergences in (|98l) . 

We recall relation (fT3"|) . We need to expand the various 
terms in the counterterm action above around D = 4, by 
making use of R = (D-l)(D-2e)H 2 and equation (fj"0"2]) : 



d 2 V 

dip 2 



R ~ 6(2-e)H 2 + (Z?-4)(7-2e)£f 2 (105a) 
(tp ) ^ 2(3 - t)eH 2 + (D - A)2eH 2 (105b) 



Kg^d^ipdvip ~ ~4,eH 2 - (D - 4)2eH 2 



(105c) 



where we neglected all 0{(D — 4) 2 } contributions in the 
equations above. We can now easily solve for the coeffi- 
cients in the counterterm action: 



Cl 



<■-> 



C3 



C4 



\d\ 2 /x 1 



48tt 2 D-A 

2 ,,D-4 



\d\ 2 n 



* Y 8tt 2 D - 4 2 
2 \c\ 2 -\d\ 2 /i p ~ 4 f 
~ 69y 32tt 2 D-i +C ' 3 
, \c\ 2 - \d\ 2 ^ , j 



-9y 



in 2 D-A 



C4 



(106a) 
(106b) 
(106c) 
(106d) 



where the divergent coefficients in D = 4 are fixed to 
cancel the divergences occurring in (|98|) . and where e.g. 
c\ is a finite but arbitrary coefficient of the counterterm 
action. The rcnormalised one loop effective action thus 
reads: 
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(N 2 -MI 2 ) 



i m 2 F 2 (l - e) 2 + I m 4 



£T 2 (l-e) 5 



2tt/i 2 



_ m 2 ff 2(7_ 5e ) 



12 



+ J dm (mH 2 {l - e) 2 + m 3 ) ^'( 2 ~ <) + V>(2 + iQ 
+2H 3 (l-e) 3 lm(cd*) J dmT(2 + iQT(2 - <)| 



(107) 



The effective action above comprises the one loop back- 
reaction of fermions. This term gives rise to additional 
physical corrections to the (classical) Friedmann equa- 
tions, thus obtaining the so-called quantum corrected 
Friedmann equations. We will study its effect on the 
evolution of our universe and its impact on the evolution 
of the scalar field to which the fermions are coupled in 



r 



a future publication. In this equation we have neglected 
all contributions to the one loop effective action arising 
from the trace anomaly. 



If we expand equation (|107[) to first order in e, we can 
find out how large the first order correction to the de 
Sitter result is. We have: 



l.rcn 



I d^X^J-\(\c\ 2 -\d\ 2 ) 



dmi m H 



m 2 H 2 \ ^ - 7e - In 



H 2 



2tt(i 2 



—m 

2 



(108) 



H 



V>'(2 + ifhJH) - 7/7(2 - im/H) 



2mH 2 



ip(2 + im/H) + ip(2 - im/H) 



2H 3 lm(cd*) / dmT{2 + im/H)T{2-im/H) 



Tfl / 

i— ip(2 + im/H) - xp(2 - im/H) 
H V 



O e 2 ) , 



where the prime denotes a derivative with respect to 
the argument, and the superscript A represent the well- 
known de Sitter result, which trivially follows from (|107p 
by letting e — > 0. 

Let us at this point however already consider the f 3> 1 
limit, to extract the physical behaviour of the effective 
potential resulting from (|107[) . Note firstly that: 



rv 



d 4 xy^(-V eS (</>)) 



(109) 



where V e s(4>) denotes the one loop effective potential for 
4>. In order to evaluate the integrals over m occurring in 
equation (|107p , we need the following asymptotic expan- 
sions ((€»): 



lim r(2 + iQT(2 - iQ = 27r|C| 3 e _7r|cl (110) 
C — >o ° 



l + T 2+ °{^ 



and also: 



lim [ii{2 + iQ+i>{2-iQ] 



ln[C 2 



13 1 

YC 2 



(111) 



60 C 4 U 6 



Keeping the leading order terms of the effective potential 
(TlTO)) in this limit only, we find from (fTU7|> and (jllOHllip : 



Kff(</>) 



\4- 



16tt 2 
Im(cd*) £ 



<?> 4 ln 



2i2 



2nfi 2 



(112) 



g^ 3 H(l-e)exp 



9y\4>\ 
H\l -< 



where equation (TlU)) has been used. Without study- 
ing the dynamics resulting from this asymptotic effective 
potential, let us briefly consider its stability properties. 
Clearly, this potential is unstable in the vacuum (|d| = 0) 
because the effective potential is unbounded from below 
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in this case. This wild behaviour has already been recog- 
nised in de Sitter spacetime [l9l ]. 

An interesting phenomenon can be observed precisely 
when | c | 2 = |rf| 2 = 1/2. This is realised for example in 
a thermal state with a temperature T much larger than 
the UV cutoff of the theory. In this case, the asymp- 
totic effective potential is exponentially suppressed and 
its (small) contributions can either be positive or nega- 
tive. 



E. Discussion of the Renormalised Effective Action 

Let us make a few general remarks. Consider quantum 
field theories in de Sitter spacetime. The Hubble param- 
eter H in de Sitter spacetime is a constant. An effective 
action for fcrmions or scalars in de Sitter spacetime con- 
tains divergences for example proportional to a power of 
H (see e.g.: [H Q]| [H, [Hj] ) . The same statement holds 
for divergences in the stress energy tensor in for example 
the trace anomaly literature (see e.g.: [2f| At this 

point, one can interpret the results differently: 

1. The Hubble parameter can quite conservatively be 
treated as a constant. The argument for this case 
is simple: throughout the calculation one simply 
assumed a constant Hubble parameter. 

2. The Hubble parameter is promoted to a dynamical 
quantity: 



H -► H(t) . 



(113) 



A reason for supporting this case would be the fol- 
lowing: when e <C 1 one would find the de Sitter 
result as a leading order contribution in each time 
interval where the Hubble parameter varies only 
adiabatically slowly, a spacetime also known as lo- 
cally de Sitter spacetime. 

Depending on which of the two interpretations one ad- 
heres to, one would renormalise differently. If one advo- 
cates the first point of view, all terms of e.g. an effective 
action merely contribute to the renormalised cosmolog- 
ical constant. If, however, one prefers the second inter- 
pretation, one adds (dynamical) curvature invariants in 
order to cancel these divergences in a covariant manner. 

We observe the following. In the effective action (11071) 
above, the Hubble parameter is a dynamical quantity. 
Hence, covariant counterterms added to cancel UV di- 
vergences can and should be treated dynamically. This 
allows us to study the backreaction of fermions in a dy- 
namical manner. 

Albeit strictly speaking not allowed to relax the as- 
sumption H = const in de Sitter spacetime, we conclude 
that promoting the Hubble parameter to a dynamical 
quantity H = H(t) has been a useful approach. When 
in our calculation we keep e constant as assumed, result 
favours the second interpretation above. 



Let us now return to our calculation to see what the 
above considerations can tell us about our effective action 
([§3"| . We face the following two options: 

1. We treat e as a constant. 

2. We promote e to a dynamical quantity, as already 
stated in (|103p . We could call this a spacetime with 
an almost constant deceleration. 

Despite the fact that throughout the calculation we 
have assumed that e is a constant, it will come as no 
surprise that we argue in favour of the latter option for 
two reasons: if one would allow e to change, one would 
probably find our main result for the effective action be- 
fore renormalisation (|93p in each time interval where e 
changes only very slowly. Secondly, the latter option is 
the generalisation of the confirmed de Sitter case above. 

An obvious disadvantage of our choice is that we can 
only renormalise on-shell. Consider the divergences in 
As already mentioned, it is not possible to identify 
local covariant counterterms of curvature invariants only 
that remove all singularities in H{t) and e(t) at the level 
of the effective action. Therefore, we relied on an on- 
shell renormalisation technique and inserted the zeroth 
order (tree level) Fricdmann equations of motion into the 
effective action at first order in h. 

Finally, let us try to comfort the reader who would 
rather keep e constant, but who would simultaneously 
like to understand how to improve on his final result for 
the renormalised effective action. In this case, the only 
appropriate counterterms to renormalise (|98|) available to 
us are the following: 



^alt 



(114) 



We can then renormalise the theory as usual 3 , with 
„ . 2\c\ 2 ~\d\ 2 M°- 4 (1-c) 2 



24tt 2 D-4 2-e 



ft 



-fh 



\d\ 2 // 



8ir 2 



D 



0t 



af (115a) 
(115b) 



We could now study mode mixing 4 . One could simply 
allow e and m/H to vary in time and choose (if neces- 
sary infinitesimally small) time intervals in which e and 



3 The case e = 2 would have to be considered separately, for R = 
in a radiation dominated universe. 

4 As already formulated in [56ll a (numerical) study of mode mixing 
is the answer to the fundamental problem that we can only solve 
for the mode functions exactly for a limited number of choices of 
a(t). In their paper, Tsamis and Woodard study mode mixing 
for a massless minimally coupled scalar field in constant e FLRW 
spacetimes. The basic strategy is as follows: in each - possibly 
infinitesimally small — time interval e can be well approximated 
by a constant. The mode functions in two adjacent time inter- 
vals can be related by imposing two conditions: (i) continuity of 
the mode functions and (ii) continuity of their first derivatives. 
Hence, the mode functions at some later time r\ can be obtained 
from the initial condition at r/ by means of a transfer matrix. 
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m/H vary only slowly. Subsequently, one can match the 
solutions in two neighbouring intervals and study mode 
mixing in the spirit of [56[. In this case, this procedure 
would imply that we would have to renormalise for each 
time interval separately, for the coefficients in equation 
(|115|) depend on e. This is not surprising: a sudden jump, 
even infinitesimally small, in for example m/H is not a 
physical phenomenon and would call for non-local coun- 
terterms (see e.g.: [57|). 



VI. CONCLUSION 



In this paper, we have constructed the fcrmionic prop- 
agator in realistic FLRW spacetimes with constant decel- 
eration parameter q = e — 1. Moreover, we assumed that 
m/H is also constant, which can be realised in Yukawa 
theory by means of a scalar field for which (j> cx H. We 
have derived the propagator both in the (trivial) massless 
case (f2"5j) and in the massive case (|78|) . 

Two pillars support our derivation. We split the 
fermionic degrees of freedom into a direct product of chi- 
rality and helicity eigenspinors. We normalise the spinors 
using a consistent canonical quantisation and we require 
orthogonality of particle and anti-particle spinors. 

Moreover, we verify that upon sending e to zero, we 
recover the known vacuum-to-vacuum de Sitter results 
in the literature. 

We calculate the one loop effective action induced by 
fermions using our propagator. The one loop backreac- 
tion arises from integrating out a free, quadratic fermion 
field and, using dimensional regularisation, this generates 
a correction to the (classical) Friedmann equations. This 
effective action allows for the first time for a dynami- 
cal interpretation of H such that the study of dynamical 
backreaction cither on the background spacetimc or on 
scalar fields coupled to these fermions is within reach. 
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Appendix A: PROPERTIES OF HANKEL 
FUNCTIONS 



In section IIV A 41 we used the following well-known 
identities involving Hankel functions: 

H { _%) = e^HP(z) (Ala) 

H^ v (z) = e-^Hl 2 \z) (Alb) 

{H^{z)Y = H$(z*) (Ale) 

HP(e*"z) = -H^l{z) = -e-^ v Hi 2 \z) (Aid) 

H^ (e~™z) = -H { H{z) = -e^HP(z) . (Ale) 
The Wronskian of two Hankel functions reads: 



4/ 



W[HP(z),hW(z)} = - 



(Alf) 



Moreover, we made use of the following recurrence rela- 
tion: 

H^{z) = ^H^{z) + -H^{z). (Alg) 

dz z 



Appendix B: FOURIER TRANSFORMING 
HANKEL FUNCTIONS 



Consider the following integral equation: 



1 = 



dec x'" +1 J tl (cx)K v {ax)K^(bx) . 



(Bl) 



Here, J^z) is a Bessel function of the first kind and 
K v {z) is a MacDonald function. Note that for all in- 
tegrals in ([76]) both a and b are purely imaginary. The 
requirements from [46| are such that: 



Re(a) > 

Re(6) > 

c > 

Re(/i±^) > -1 

Re(/i) > -1 



(B2a) 
(B2b) 
(B2c) 
(B2d) 
(B2e) 



Hence we introduce small real contributions in the argu- 
ments of both MacDonald functions only to take these to 
zero again at the end of the calculation. We thus write: 



a = 5 
b = 5 



ict 
i(3. 



(B3a) 
(B3b) 



We choose 5 > to make the integral convergent. We 
now use equation (6.578.10) from [46j to find: 



_ ^r( M + ^ + i)r(/i-^ + i) „_ 



where: 



6 2 + c 2 



2ab 



'(«). (B4) 



(B5) 
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Here, P _i 2 (u) represents the associated Legendre func- 
v . 2 

tion which we can rewrite in terms of a Gauss' hyperge- 
ometric function %Fi by means of equations (8.702) and 
(9.131.1) of 0: 

4(o6)"+ 1 r0i + 3/2) 1 ' 

( 3 1-u 

x 2 Fi \fj, + i/ + l,fj, — v + l;fj,+ -; — - — 

Upon taking 8 — ► we find: 
1 1 



(B7) 



lira 

5 



Finally note that we can further simplify the argument 
of the hypergeometric function in equation (|B6p to: 



c 2 - (a + (3 + %8f 
4a/? 



(B8) 



Let us finally stress that we have introduced 8 > to 
make the integral convergent. It is however of significant 
physical importance for it is the same 8 that dictates the 
e pole prescription in the propagator (see equations ([751 - 
[79[) and ([22)) ). This should not come as a surprise: we 
observe the same behaviour in e.g. Minkowski and de 
Sitter spacetimes. 
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